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We  consider  the  instability  of  Couette  flow  between  rotating  concentric 

cylinders  of  infinite  length.  Let  R  ,  R^  and  flj,  denote  the  radii  and 

angular  velocities  of  the  inner  and  outer  cylinders,  respectively;  and  let 

I,  2 

n  “  R^R.,  and  u  =  Hg/fij.  We  also  introduce  the  Taylor  number  T  *  -UAfl^d  /v 
with  a  -  (n2 r22  -  o1r18)/(r22  -  Rj2). 

It  is  knovn  from  the  linear  stability  analyis  of  Krueger,  Cross,  and  DiPrima 

(1966)  for  n  near  one  that  (i)  for  11  s  0  and  increasing  T,  Couette  flow' 

first  becomes  unstable  to  an  axisymmetric  disturbance  at  I  '  Tc(m  *  0,  n), 

(ii)  in  the  absence  of  an  axisymmetric  disturbance  Couette  flow  would  become 

unstable  to  a  non-axisymmetric  disturbance  at  a  value  of  T  just  slightly 

greater  than  T  ,  and-(iii)  at  u  sufficiently  negative  the  critical  value  of 
c 

T  occurs  for  a  non-axi symmetric  disturbance.  Thus,  for  fixed  n  near  one  there 
are  parameter  values  (w^,  T^)  for  which  axisymmetric  and  non-axi symmetric 
disturbances  are  simultaneously  critical  for  Couette  flow.  We  wish  to  study 
the  bifurcation  problem  in  the  neighborhood  of  (t>1,  T( ) ,  where  we  note  that 
instability  to  an  axisymmetric  disturbance  gives  a  steady  bifurcation  while 
instability  to  a  non-axisymmetric  disturbance  gives  a  Hopf  bifurcation. 

In  order  to  develop  a  mathematically  tractable  problem,  we  assume  that  the 
axial  wavenumber  of  the  axisymmetric  and  non-axisymmetric  disturbances  are 
the  same,  and  we  choose  the  wavenumber  to  be  that  corresponding  to  critical 
conditions  for  an  axisymmetric  disturbance.  We  let  X  denote  the  dimensionless 
axial  wavenumber  (scaled  with  respect  to  Rj-R^ )  and  m  denote  the  number  of 
waves  in  the  azimuthal  direction.  Or.  Peter  Eagles  has  earried  out  the 
necessary  calculations  for  n  “  0.9$  to  determine  the  points  (w,,  T1 )  for 
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simultaneous  instability  to  an  axisymmetric  disturbance  (\,0)  and  a  non-axi- 
syrmetric  disturbance  (X,m)  vith  the  results  given  in  the  Table. 


1 

2 

U 


1 

-0.73976 

-0.75123 

-O.8028I4 


Table 

T1 

11.973.P7 

12,319.01 

1»i,015.01 


X 

3.  *182 
3-507 
3-628 


Thus,  for  example,  for  n  =  O.95  and  u  ■  -0.73976,  Couette  flow  is  simultaneous¬ 
ly  unstable  to  an  axisymaetric  disturbance  vith  wavenumbers  (3. *>82,0)  and  a 
non-axi symmetric  disturbance  with  wavenumbers  (3.1<82,1 ).  We  alsonotethat  at 
(v^ ,  T1)  there  are  6  critical  modes  vith  axial  (Z)  and  azimuthal  (6)  depend¬ 
ence  as  follows:  cos  XZ,  sin  XZ,  exp  (±  imB)  cos  XZ  and  exp  (i  imB)  sin  XZ. 

The  first  step  in  the  bifurcation  analysis  is  to  represent  the  components  of 
the  disturbance  by 


A  u  (r)  cosXZ  +  A  u  (r)  sinXZ  + 
c  c  s  s 

Bcvc(r)  exp  (im8)  cosXZ  +  B^v^r)  exP  (£■•)  sinXZ  +  (l) 

S^^tr)  exp  (-im0)  cosXZ  ♦  S^v^ ( r )  exp  (-imB)  sinXZ  +  T 

Here  the  functions  uc,  u^ ,  v^ ,  denote  the  radial  dependence  of  the  res¬ 
pective  modes.  The  quantities  A  ,  A  ,  B  ,  B  are  scalar  functions  of  time 

c  s  c  s 

vith  Ac>  Ag  varying  in  the  real  numbers  and  Bc>  B&  in  the  complex  domain. 
The  residual  component  T  is  orthogonal  to  the  six  critical  modes. 

We  are  now  able  to  reduce  the  nonlinear  partial  differential  equation  to  a 
finite-  dimensional  center  manifold  of  dimension  six;  this  reduction  is 
performed  by  expressing  T  in  terms  of  the  coefficients  A^ ,  A# ,  Bc>  B^.  The 
nonlinear  ordinary  differential  equations  for  these  coefficients  describe 
the  'flow*  in  the  center  .manifold  and  are  similar  to  those  derived  by  Davey, 
DiPrima  and  Stuart  (1968)  for  the  case  v  “  O. 

We  give  particular  attention  to  the  interaction  of  the  three  modes  cosXZ 
and  exp  (ZimB)  sinXZ,  vhich  can  give  rise  to  Taylor-vortex  and  vavy-vortex 
flows  vhich  are  often  observed  in  experiments.  So  ve  put  A^  *  B^  *  0  and 
we  obtain  the  following  equations  where  terms  of  the  order  7  have  been 
omitted: 
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dAc 

dt 

dBs 

dt 


qAc  *  aJAc3  +  a^AjBj2  ♦  a^5  +  *gAc3lBs|2  +  a^jBj*. 


V.  +  b,BJBsl2  +  Wc* +  v.V +  Vs>B,l: V  + 


(2) 


Complex  conjugate  equation  for  §s- 

The  scalar  quantities  a.  are  real  and  the  b.  are  complex  numbers.  In  the 

J  J 

sequel  ve  shall  denote  by  b.  and  b..  the  real  and  imaginary  parts  of  b.. 

Jr  jx  J  * 

After  performing  the  transformation  Ac  “  x,  B#  =  y  exp  iCb^t+O  the  equation 
for  ♦  uncouples  from  those  for  x  and  y  and  ve  obtain 


dx  3  2  5  3  2  *• 

dt  "  V  -  V  *  V*  +  V  “ 

dt 


+  agx  y  +  a^xy 
b0r*  +  bir*3  +  b4rx2y  *  V**  *  +  b7ry? 


(3) 


The  coefficients  a^  and  b-r  depend  on  V  and  T,  vhile  aQ(|ij,  T1 )  *  bQr(u1 ,  T1) 

=  0. 

It  is  now  possible  to  study  changes  in  the  phase  plane  as  a  function  of  the 
tvo  parameters  ir-tij  and  T-T^ .  It  is  natural  first  to  analyze  system  (3)  vith 
the  fifth  order  terms  left  out.  Keener  (1976)  and  Langford  and  Iooss  (1980) 
consider  pairs  of  coupled  ordinary  differential  equations  vith  the  same  cubic 
nonlinearities  as  (3)  by  taking  a  ,  a^,  blr,  b|jr  fixed  and  varying  aQ  and  bQr . 
The  nature  of  the  changes  in  the  phase  plane  depends  of  course  on  the  specific 
values  of  the  1  order  coefficients  and  the  cubic  systems  of  o.d.e.'s  can 
be  classified  accordingly. 

Vhile  the  investigations  by  Keener  and  Langford-Iooss  vere  carried  out  for 
arbitrary  nonvanisbing  values  of  the  3rd  order  coefficients,  it  turns  out 
that  in  tbe  problem  vhich  is  presently  under  discussion  some  of  these 
coefficients  are  very  small  and  change  sign  for  values  of  w,  T  very  close 
to  |i.j,  Tj.  For  example,  if  ve  let  v  take  values  dose  to  y^  and  -compute 
T  ■  T  (v,n  “  -95,  1  “  3.M52,  n  «  0)  then  a,(y,T)  -  0  for  v  *  . 68939 •  vhich 

is  indeed  very  close  to  u1  *  -.73976.  For  this  reason  tbe  analysis  of  (3) 
has  to  be  taken  to  fifth  order  and  the  dependence  of  a^ ,  a^ ,  blr»  on  v 
and  T  has  to  be  taken  into' account  up  to  first  order  in  y-y^  and  T-T^. 

Tbe  numerical  evaluation  of  tbe  very  complicated  expressions  for  tbe  fifth 
order  coefficients  have  been  carried  out  by  Dr.  Eagles  and  tbe  analysis  of 
tbe  changes  in  tbe  phase  plane  is  in  progress. 
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The  fall  system  of  six  amplitude  equations  is  also  studied,  tout  only  at 
cubic  order  in  tfae  amplitudes.  For  this  case  there  does  not  appear  to  toe  a 
transformation  that  vill  reduce  the  order  of  tbe  system. 

For  the  case  u  *=  0,  ve  have  noted  that  the  axisymmetric  and  non-axisymaetric 
modes  are  not  simultaneously  unstable.  However,  it  is  possible  to  imbed  this 
problem  in  a'larger  class  of  mathematical  problems  which  ve  can  analyze 
using  the  techniques  of  center  manifold  theory.  In  this  way  ve  present  a 
rational  basis  (though  admittedly  complicated)  for  the  analysis  of  the 
successive  bifurcations  from  Couette  flow  to  Taylor-vortex  flow  and  then  to 


vavy-vortex  flow  for  increasing  T  with  y  «  0  as  first  studied  toy  Davey, 
DiPrima,  and  Stuart  (1966). 
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